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Topological crystalline insulators (TCI) are new topological phases of matter protected by crystal 
symmetry of solids. Recently, the first realization of TCI has been predicted and observed in 
IV- VI semiconductor SnTe and related alloys Pbi_^Sn I (Te, Se). By combining k ■ p theory and 
band structure calculation, we present a unified approach to study topological surface states on 
various crystal surfaces of TCI in IV- VI semiconductors. We explicitly derive k ■ p Hamiltonian 
for topological surface states from electronic structure of the bulk, thereby providing a microscopic 
understanding of bulk-boundary correspondence in TCI. Depending on the surface orientation, we 
find two types of surface states with qualitatively different properties. In particular, we predict that 
(111) surface states consist of four Dirac cones centered at time-reversal-invariant momenta f and 
M, while (110) surface states consist of Dirac cones at non-time-reversal-invariant momenta, similar 
to (001). Moreover, both (001) and (110) surface states exhibit a Lifshitz transition as a function 
of Fermi energy, which is accompanied by a Van-Hove singularity in density of states arising from 
saddle points in the band structure. 
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Structure and symmetry play an important role in 
shaping electronic properties of periodic solids. It is a 
common phenomenon that materials made of chemically 
similar elements arranged in the same crystal structure, 
like diamond and silicon, often have qualitatively simi- 
lar electronic properties. This empirical relationship be- 
tween structure and property arises from the fact that 
the essential electronic properties of many solids are un- 
derstandable in terms of orbitals and bonds, the char- 
acteristics of which depend mostly on crystal structure. 
For example, both diamond and silicon possess sp 3 hy- 
bridized orbitals in tetragonal structure. On the other 
hand, the quantum theory of solids is based on itinerant 
Bloch waves that form energy bands in momentum space. 
The global structure of band theory allows for unconven- 
tional energy bands that are topologically nontrivial as 
a whole entity. Such band structures are fundamentally 
different from a Slater determinant of atomic orbitals, 
and give rise to topological states of matter exhibiting 
universal and quantized properties that are absent in or- 
dinary solids. Celebrated examples of such topological 
states include quantum Hall insulators [T] and topologi- 
cal insulators [2H3]. 

The topological aspect of band theory gives birth to 
interesting exceptions to the empirical rule of structure- 
property relation. We demonstrated [5] the proof of prin- 
ciple that there exist distinct classes of band structures 
within the same crystal structure, which cannot be adia- 
batically connected under deformations preserving cer- 
tain point group symmetries. Those nontrivial band 
structures are characterized by topological indices, and 
thereby define a new topological state of matter dubbed 
topological crystalline insulators (TCI). A hallmark of 
TCI is the existence of gapless surface states on those sur- 
faces that preserve the underlying crystal symmetry [5J. 



Recently, Hsieh et al predicted [7] the first material 
realization of TCI in IV- VI semiconductors SnTe and 
related alloys Pb a; Sni_ a ;(Te, Se). The nontrivial band 
topology in these materials is characterized by an integer 
topological invariant known as mirror Chern number [5], 
arising from the mirror symmetry with respect to the 
(110) plane of the rocksalt structure (and its symmetry- 
related ones). A consequence of this electronic topol- 
ogy is the existence of topological surface states on crys- 
tal surfaces that preserve at least one such mirror sym- 
metry, such as (001), (111), (HO). In particular, the 
(001) surface states consist of four Dirac cones at low en- 
ergy, which are located at generic point A on the mirror- 
symmetric line YX. Importantly, these four Dirac cones 
are spin-momentum locked with the same chirality, which 
is a unique hallmark of the TCI phase in IV- VI semi- 
conductors. These predicted surface states were subse- 
quently observed in angle-resolved photoemission spec- 
troscopy experiments on SnTe[3], Pbi-^Sn^SefTU] and 
Pb^Sni-zTepJ]. In particular, the spin-texture of these 
surface states observed by Xu et a/QT] (see also Ref.[T^]) 
provides a direct spectroscopic measurement of the mir- 
ror Chern number [T51ll4) . 

The materialization of TCI opens up a new venue 
for topological states of matter in a much larger num- 
ber of material classes than previously thought [ 15141 9| . 
thereby triggering intensive activities [20 24 . From a ma- 
terial viewpoint, IV- VI semiconductors have high mobil- 
ities and exhibit a wide range of electronic properties 
(e.g., magnetism, ferroelectricity and superconductivity) 
that can be easily tuned by alloying, doping and strain. 
The technology for synthesizing and engineering these 
materials, in both bulk and low-dimensional form, has 
been well-developed by decades of efforts. Therefore TCI 
in IV- VI semiconductors provide an extremely versatile 



platform for exploring topological quantum phenomena 
and novel device applications. 

In this work, we combine k ■ p theory and band struc- 
ture calculation to study topological surface states on 
various crystal surfaces of TCI in IV- VI semiconductors. 
We present a unified approach to derive k ■ p Hamilto- 
nian for surface states from the electronic structure of 
the bulk, thereby providing a microscopic understanding 
of bulk-boundary correspondence in TCI. We find that 
low-energy properties of surface states are solely deter- 
mined by the surface orientation, and can be classified 
into two types. In particular, we predict that (111) sur- 
face states consist of four Dirac cones centered at time- 
reversal-invariant momenta T and M|25j. while (110) sur- 
face states consist of Dirac cones at non-time-reversal- 
invariant momenta, similar to (001). Moreover, both 
(001) and (110) surface states exhibit a Lifshitz transition 
as a function of Fermi energy, which is accompanied by 
a Van-Hove singularity in density of states arising from 
saddle points in the band structure. Our results provide 
a much-needed basis for further investigations on TCI 
surface states. 

We begin by reviewing k ■ p theory for the bulk band 
structure of TCI, from which surface states are derived. 
The band gap of IV- VI semiconductors is located at four 
L points. For the ionic insulator PbTe, the Bloch state 
of the valence band at L is derived from the p-orbitals of 
the anion Te, and that of the conduction band from the 
cation Pb. In contrast, SnTe has an inherently inverted 
band ordering, in which the valence band is derived from 
the cation Sn and the conduction band from Te. This 
band inversion relative to a trivial ionic insulator gives 
rise to the TCI phase in SnTe [7J. 

The band structure near each L point can be described 
by k ■ p theory in the basis of the four Bloch states at L, 
ifccr 8 (L), where a = 1(— 1) refers to the state derived from 
the cation (anion), and s labels the Kramers degeneracy. 
The k ■ p Hamiltonian H(k) is given by (see [7J and ref- 
erences therein): 



.ff(k) = ma z + va x (k 1 s 2 - k 2 s\) + v'k 3 a y 



(1) 



where k 3 is along the TL direction, and k\ is along the 
(llO) axis of reflection, a and s are two sets of Pauli 
matrices. The sign of m in (fTl) captures the two types of 
band ordering: in our convention m > for PbTe and 
m < for SnTe. 

Electronic structures of TCI surface states depend cru- 
cially on crystal face orientations. We find surface states 
with qualitatively different properties on two types of 
crystal surfaces. For type-I surface, all four L points are 
projected to different time-reversal-invariant momenta in 
the surface Brillouin zone. This is the case for (111) sur- 
face, for which L\ is projected to T and (L 2 ,L3, L4) are 
projected to three M points. For type-II surface, differ- 
ent L points are projected to the same surface momenta. 
This is the case for (001) surface for which (Li, L 2 ) —> X\ 




FIG. 1: Bulk Brillouin zone of rocksalt structure and its pro- 
jection to (001), (111) and (110) surface Brillouin zone. For 
(001) surface, L\ and L2 are both projected to X\, and L3 
and Li are projected to X2; For (111) surface, L\ is projected 
to F and the other three L points are projected to M points; 
For (110) direction, L\ and Li are projected to X and L3 and 
Li are projected to R. The shaded plane passing through T, 
L\ and L2 in 3D Brillouin zone is invariant under reflection 
about (110) plane in real space. 



and {Ls,Li) — > X 2 , as well as (110) surface for which 
(Li,L 2 ) —> X and (£3,1,4) — > R. The projection from 
bulk to surface Brillouin zone is illustrated in Fig.l. 

(Ill) surface: type-I surface states can be obtained 
straight-forwardly from the continuum k ■ p Hamiltonian 
Q. Following the spirit of Ref.[2B], we model vacuum 
as a trivial insulator (like PbTe) with an infinite gap 
m = M > 0, where M —> +00. Surface states can now 
be obtained by solving a domain wall problem in which 
the Dirac mass m changes sign across the interface be- 
tween TCI and vacuum. It is well-known from field the- 
ory that two-dimensional massless Dirac fermions form at 
such an interface [2"TH2l?] . Due to the presence of four L- 
valleys, surface states consist of four copies of such Dirac 
fermions. These four Dirac fermions are located at four 
distinct 2D momenta that correspond to the projection 
of the four L points onto the type-I surface. Because 
of their different in-plane momenta, the four L-valleys 
cannot couple with each other as long as in-plane trans- 
lation symmetry is present, and therefore independently 
give birth to four branches of Dirac surface states. 

A prime example of type-I surface is (111). It follows 
from the above analysis that (111) surface states consist 
of four Dirac cones: one at f , and three others at M. 
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FIG. 2: Band structure of (111) surface for Sn and Te termi- 
nation. Four Dirac pockets are present: one at F and three at 
M. This leads to two counter-propagating states with oppo- 
site mirror eigenvalues on TM, as predicted from the mirror 
Chern number Um = —2 in the TCI phase]?]. 



The k ■ p Hamiltonians at T and M are given by 



.fff(k) = v{k\s 2 — k 2 s\) 
H Ai(k) = v 1 k 1 s 2 - v 2 k 2 s 1 



(2) 



where k\ is along TK direction, and k 2 is along the 
mirror-invariant TM direction. The presence of these 
Dirac pockets is confirmed by our band structure calcu- 
lations on SnTe, based on the tight-binding model [50], 
see Fig. 2. The Dirac points are found to lie close to 
the top (bottom) of the valence (conduction) band for 
Sn (Te) termination. Such surface states are qualita- 
tively similar to interface states between PbTe and SnTe 
studied in early theoretical works [2TH2T)] . The advent of 
TCI has revealed that such interface states are topologi- 
cally equivalent to surface states of SnTe (but not PbTe), 
and their robust existence is topologically protected by 
the (110) mirror symmetry. This protection can be un- 
derstood from the two branches of counter-propagating 
surface states on the mirror-symmetric line TM. Our 
tight-binding calculation shows that these two branches 
have opposite mirror eigenvalues and therefore cannot 
couple with each other to open up a gap, as long as mirror 
symmetry is preserved. This unusual property of (111) 
surface states agrees with the prediction based on mir- 
ror Chern number]?], and constitutes another important 
hallmark of TCI. Since mirror eigenvalue of a spin-1/2 
fermion is intimately related to its spin state, it follows 
that spin-textures on all four Dirac pockets, at both T 
and M, have the same chirality. 

(001) surface: type-II surfaces are much more in- 
teresting. In this case, two L points are projected to 
the same momentum in the surface Brillouin zone. As a 
result, they interact with each other to create topologi- 
cal surface states that are different from type-I surfaces. 
This interaction between L- valleys arises from physics at 
the lattice scale, which is not captured by the previous 
field-theoretic approach |27H29j . For example, the previ- 
ous ab-initio calculation 7] shows that (001) surface has 
a pair of Dirac cones away from X (the projection of L), 
which cannot be understood from field theory alone. 

To obtain type-II surface states, we proceed in two 



steps. First, we study a hypothetical smooth interface 
between TCI and a trivial insulator, in which the Dirac 
mass is slowly-varying in space and gradually changes 
from negative to positive over many lattice constants 
across the interface. Scattering between the two L val- 
leys projected to the same momentum on the surface re- 
quires large-momentum transfer in the direction normal 
to the interface and hence vanishes in this smooth limit. 
We are thus justified to treat different valleys indepen- 
dently and use the continuum field theory to solve for 
interface states. The effective Hamiltonian that we de- 
rive for these interface states serves as a starting point for 
surface states. Next, imagine deforming the smooth in- 
terface into the atomically sharp surface, which adiabati- 
cally connects interface states to the actual surface states. 
This deformation procedure introduces inter-valley scat- 
tering processes at the lattice scale, which are represented 
by additional terms in the effective Hamiltonian for the 
surface. Such terms must satisfy all the crystal symme- 
try at the surface, and therefore can be enumerated by 
a symmetry analysis. By incorporating these terms into 
the surface Hamiltonian derived in the previous step, we 
obtain the final form of k • p theory for type-II surfaces. 
We now apply this approach to study (001) surface. 
Starting from a smooth interface, there exist two mass- 
less Dirac fermions at Xi arising from the L\ and L 2 
valley respectively, and likewise for the symmetry-related 
point X 2 . These two flavors of Dirac fermions have iden- 
tical energy-momentum dispersions, resulting in a two- 
fold degeneracy at every k. The k ■ p Hamiltonian for 
this smooth interface is given by 



"X\ \ ) — \ v xk x Sy 



VyfVySx) 



(3) 



where the momentum {k x , k y ) is measured from X\, with 
k x parallel to TX 2 and k y parallel to TXi . Here I is iden- 
tity operator in the flavor space, and s is a set of Pauli 
matrices associated with the two components of each fla- 
vor that are intimately related to electron's spin. The 
velocities in x and y directions are generically different. 
Next, we perform a symmetry analysis to deduce the 
form of those additional terms associated with physics at 
the lattice scale, which need to be added to ([3]). Note 
that X is invariant under three point group operations: 
(i) x — » ~x reflection (M x ); (ii) y — > —y reflection (M y ); 
(iii) two- fold rotation around surface normal (C 2 ). These 
symmetry operations, plus time reversal transformation 
0, are represented by the following unitary operators in 
our k ■ p theory: 



M x 

My 

C 2 

e 






ISyK 



(4) 



Here M x preserves the L\ and L 2 valley in the bulk and 
only acts on electron's spin, whereas both M y and C 2 in- 



terchange L\ and L2 and hence involve a flavor-changing 
operator r x . To zero-th order in k, we find two symmetry- 
allowed operators: t x and t v s x . Therefore our k-p Hamil- 
tonian for (001) surface states is given by: 



H Xi (k) = {v x k x s y ~ Vyk y s x ) + mT x + Ss x T y 



(5) 



Note that the two additional terms, parameterized by to 
and 5, are off-diagonal in flavor space, which correctly 



describe inter-valley scattering at the lattice scale. 

The k ■ p Hamiltonian ([5]) is a main result of this work. 
We now show that Hx (k) captures all the essential fea- 
tures of type-II surface states. By diagonalizing H Xl (k), 
we obtain four surface bands with energy-momentum 
dispersions Eh (k), — £#(k), £x(k) and — -Ex(k) respec- 
tively, where Eh,lO*) is given by 



E 



H,L 



(k) = Jm 2 + 5 2 + v' 2 k x + v 2 k 2 ± 2Jm 2 v 2 k x + (to 2 + S 2 )k 2 v, 



(6) 
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FIG. 3: k-p band structure for (001) surface states. A pair of 
low-energy Dirac cones, located at Ai and A2 on the line XT, 
is formed by the interaction between two high-energy Dirac 
bands centered at X. k-p parameters are obtained by fit- 
ting with ab-initio results (shown by red dots) on SnTe taken 
from Ref.[7]: v x = 2.4eV-A, v y = 1.3eV-A, m = 70meV, and 
8 — 26meV. Constant-energy contour evolves rapidly with 
increasing energy from the Dirac point, changing from two 
disconnected electron pockets to a large electron pocket and 
a small hole pocket via a Lifshitz transition. At this transition 
point, a saddle point S on the line XM leads to a Van-Hove 
singularity in density of states at energy E = 8. 



The resulting surface band structure is plotted in Fig. 3. 
We find two high-energy bands ±Eh which start from 
energy E x = \/m 2 + 6 2 at X and coexist in energy with 
bulk bands, as well as two low-energy bands ±El which 
mostly lie within the band gap. 

It is clear from (f5l [6]) that the two terms m and 6 arising 
from the lattice scale play a key role in shaping the band 
structure in Fig. 3. First, a finite to turns the two flavors 
of massless Dirac fermions into a "bonding" (t x = 1) 
and an "anti-bonding" (t x — —1) Dirac cones, which are 



centered at X and have energy ±to respectively. In the 
absence of 8, the lower band of the upper Dirac cone 
and the upper band of the lower Dirac cone cross each 
other at E = over an elliptical contour C in momentum 
space defined by v x k 2 +v 2 k 2 = to 2 . Next, a small <5 turns 
this band crossing into an anti-crossing via hybridization. 
Importantly, the hybridization matrix element depends 
on the direction of k, and leads to a p-wave hybridiza- 
tion gap over the contour: A(k) = 2<5 • v x k x /m, k G C. 
Importantly, A(k) vanishes along the mirror-symmetric 
line TXi corresponding to k x = 0. This is a consequence 
of the unique electronic topology of the TCI protected 
by mirror symmetry. As can be seen from Q and (pi), 
the two low-energy bands ±El have opposite M x mirror 
eigenvalues on the k y line XiT, but identical M y mirror 
eigenvalues on the k x line X\M, As a result, hybridiza- 
tion is strictly forbidden on A]T, but allowed on X\M. 
The presence of such a protected band crossing on X\T , 
but not elsewhere, leads to a pair of zero-energy Dirac 
points A12 located symmetrically away from X\ at mo- 
menta A 1.2 = (0, ±-\/to 2 + S 2 /v y ). By linearizing band 
structure near each A, we obtain the two-component 
massless Dirac fermion of Ref . [7] 



H K (6k) = v x Sk x a y 



v y Sk y cr x 



(7) 



where £k = k — A and the Dirac velocity along f X\ is 
reduced from v x : v x = v x 8/\/m 2 + S 2 . 

Our k-p theory also captures essential features of (001) 
surface states at higher energy, found in the previous ab- 
initio calculation [7]. Fig. 3 shows the band dispersion and 
constant energy contours of the k ■ p model. The sur- 
face band structure evolves rapidly with increasing en- 
ergy away from the Dirac point. For E < 6, the Fermi 
surface consists of two disconnected small Dirac pockets 
outside X. As E increases, the two electron pockets first 
grow in size, and then touch each other to transform into 
a large electron pocket and a small hole pocket, both 
centered at X. This change of Fermi surface topology 
(Lifshitz transition) occurs at a saddle point in the band 
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FIG. 4: Band structure of SnTe (110) surface, from our ab- 
initio calculations. The inset shows the surface Brillouin zone. 
A pair of Dirac cones is present on the line X — T — X, but 
absent along other high-symmetry lines. 



structure located at momentum S = (m/v x ,0) and en- 
ergy \E S \ — 5. The band dispersion near the saddle point 
is given by 



E(Sk) = E- 



2m r 



9m 



(8) 



where <5k = k — S and the effective mass tensor is de- 
fined by m xx — S/v x and m vy — m 2 /(S • v 2 )- This sad- 
dle point results in a Van-Hove singularity in density of 
states, shown in Fig. 3. 

We also compare the band structure from k ■ p the- 
ory ^ with the ab-initio result on SnTe (001) surface[7], 
shown by the fitting in Fig. 3. The quantitative agree- 
ment clearly demonstrates the validity of our k ■ p theory 
in a wide energy range measured from the Dirac point. 
We point out that the fitting can be further improved 
by further including k-linear inter-valley terms in k ■ p 
theory, in addition to the zero-th terms considered so far 
(see appendix). This leads to a more sophisticated k- 
Hamiltonian with seven independent parameters, which 
is related to a recent study by Fang et al 35j . We note 



that these k-linear terms merely modifies the result of 
H x (k) quantitatively, but do not generate the double- 
cone band structure on their own. 

(110) surface: we end by briefly discussing another 
type-II surface, (110). In this case, L\ and L 2 are pro- 
jected to X, and L 3 and L A are projected to R. Bulk- 
boundary correspondence based on the electronic topol- 
ogy of TCI predicts the existence of a pair of counter- 
propagating states with opposite mirror eigenvalues on 
rX[7]. This is confirmed by our ab-initio calculation [3Tj 
for SnTe (110) in Fig. 4, showing a pair of Dirac cones on 
the line X — T — X, but not along other high-symmetry 
lines. Since X on the (110) surface has the same symme- 
try as X\ on (001) surface (both have two mirror planes 
plus a two- fold axis) , our k ■ p theory Q and (pi) applies 
equally well to (110) surface band structure near X. 

Note: Near the completion of this work, we be- 
came aware of a preprint on tight-binding calculation of 
(Pb,Sn)Te surface states^. We also learned of a related 
work by Y. J. Wang et al. 
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APPENDIX 

In the main text, we considered lattice-scale terms that 
are zero-th order in k in studying (001) surface states. 
For completeness, we now further include k-linear terms 
that are invariant symmetry transformations Q and ex- 
amine their effects on the (001) surface band structure. 
We find four such terms: k x s z T z , k x s y T x , k y s x r x and 
k y Ty. This leads to a more sophisticated k ■ p Hamilto- 



nian for (001) surface: 



ilv (k) — (v x k x Sy Vyk y s x ) + v x k x s y T x + v x k x s z T z v y k y s x T x + v y k y T y + ttit x + os a 

I 



(9) 



We further simplify this Hamiltonian by performing a 
unitary transformation U(9) = e ies ^ T -/ 2 , Importantly, 
the form of all symmetry operators Q is preserved af- 
ter this transformation. By choosing 6 properly, we can 



eliminate the term v x without generating any new term. 
Therefore the full k-p theory for (001) surface, up to first 
order in k, is given by 



Hx ± (k) = (v x k x s y 



VyfcyS : 



:) + (<k x 



VykySx) 



VykyTy 



Ss x 



(10) 



The first two terms in ( 10 ) are inherited from the contin- 
uum theory Q without interaction between valleys. Re- 
markably, all remaining terms are off-diagonal in valley 
space, in accordance with their origin from inter-valley 
mixing at the lattice scale. The natural separation be- 
tween the continuum and the lattice-scale is the very ba- 
sis of our theory of topological surface states in the TCI 
phase of IV- VI semiconductors. 

A recent study 35 obtained a k ■ p model for (001) 
surface states, which is equivalent to ^ after a unitary 
transformation and hence can be further simplified into 



( 10 1 by eliminating a redundant parameter. Moreover, 
the model in Ref. [35] is constructed entirely by symmetry 
analysis, and therefore the physical meaning of the basis 
states in the model is not understood. 

The k-linear term (v' x k x s v — v 'k y s x )T x -\-v'J.k y T y in the 
full k ■ p theory, which is neglected in the main text, 
plays a secondary role compared to the zero-th order term 
mT x + 8s x T y . Its main effect is to break the particle-hole 
symmetry of the surface band structure (pi [6]), without 
changing any of the essential physics. Since narrow-gap 
IV- VI semiconductors have nearly mirror-like conduction 
and valence band at low energy, it is a good approxima- 
tion to neglect this k-linear term when addressing the 
main points. We do expect that including this term will 
further improve the agreement between theory and ex- 
periment. 
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